High-fT Precession modes: Axially symmetric limit of wobbling motion 
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The rotational band built on the high- AT multi-quasiparticle state can be interpreted as a multi- 
phonon band of the precession mode, which represents the precessional rotation about the axis 
perpendicular to the direction of the intrinsic angular momentum. By using the axially symmet- 
ric limit of the random-phase-approximation (RPA) formalism developed for the nuclear wobbling 
motion, we study the properties of the precession modes in 178 W; the excitation energies, B(E2) 
and B(M1) values. We show that the excitations of such a specific type of rotation can be well 
described by the RPA formalism, which gives a new insight to understand the wobbling motion in 
the triaxial superdeformed nuclei from a microscopic view point. 

PACS numbers: 21.10.Re, 21.60. Jz, 23.20.Lv, 27.70. +q 



I. INTRODUCTION 

Rotation is one of typical collective motions in atomic 
nuclei. It manifests itself as a rotational band, a sequence 
of states connected by strong electromagnetic (e.g. E2)- 
transitions. Most of the rotational bands observed so 
far are based on the uniform rotation about an axis per- 
pendicular to the symmetry axis of axially symmetric 
deformation. The well known ground state rotational 
bands and the superdeformed rotational bands with axis 
ratio about 2:1 are typical examples of this type of rota- 
tional motion. Quite recently, exotic rotational motions, 
in contrast to the normal ones mentioned above, have 
been issues under discussion, which are generally non- 
uniform nor rotating about one of three principal axes 
of deformation, and clearly indicate possible existence 
of three-dimensional rotations in atomic nuclei. The re- 
cently observed wobbling rotational bands [ E BJl Jil 
and the chiral rotation/ vibration bands [f| lalSl are 
such examples. 

Such exotic rotations are very interesting because they 
give hints to a fundamental question: How does an 
atomic nucleus rotate as a three-dimensional object? 
They may also shed light on collective motions in nu- 
clei with triaxial deformation, which are characteristic 
in these rotational bands, and are very scarce near the 
ground state region. Although the triaxial deformation 
is crucial for those exotic rotations, it is not a necessary 
condition for three-dimensional rotations to occur. For 
example, the chiral rotation is a kind of "magnetic ro- 
tation" or "tilted axis rotation" [ljj, where the axis of 
rotation is neither along a principal axis of deformation 
nor in the plane of two principal axes, but is pointing 
inside a triangle composed of three principal axes. In 
the case of the typical magnetic rotation observed in the 
Pb region, the so-called "shears band" 01 > the deforma- 
tion is axially symmetric and of weakly oblate. Similarly, 
one can think of an axially symmetric limit of the wob- 
bling motion; the so-called "precession band" , which is 
nothing else but a rotational band excited on a high- if 



isomeric state, in analogy to the classical motion of the 
symmetric top. The main purpose of the present paper 
is to investigate the precession band from a microscopic 
view point. 

In recent publications 0, 0] , we have studied the nu- 
clear wobbling motions associated with the triaxial su- 
perdeformed (TSD) bands in Lu and Hf isotopes on the 
basis of the microscopic framework; the cranked mean- 
field and the rando m p hase approximation (RPA) 0, 
13 13 13 E 0. It has been found that RPA 
eigen-modes, which can be interpreted as the wobbling 
motions, appear naturally if appropriate mean-field pa- 
rameters are chosen. The deformation of the mean-field is 
large (ea > 0.35) with a positive triaxial shape (7 ss +20° 
in the Lund convention), i.e., mainly rotating about the 
shortest axis, and the static pairing is small {A HiP < 0.6 
MeV), both of which are in accordance with the poten- 
tial energy surface calculation [20( . It should be stressed 
that the solution of the RPA eigen-value is uniquely de- 
termined, once the mean-field is fixed, as long as the 
"minimal coupling" residual interaction is adopted (see 
Sec. 1111(1 . Therefore, it is highly non-trivial that we could 
obtain wobbling-like RPA solutions at correct excitation 
energies. However, the detailed rotational frequency de- 
pendence of the observed excitation energy in Lu iso- 
topes, monotonically decreasing with frequency, could 
not be reproduced, and the out-of-band B(E2) values 
from the wobbling band were considerably underesti- 
mated in our RPA calculation. 

By restricting to the axially symmetric deformation 
with an uniform rotation about a principal axis, the an- 
gular momentum of high-spin states is built up either by a 
collective rotation, i.e., the rotation axis is perpendicular 
to the symmetry axis, or by alignments of single-particle 
angular momenta, i.e., the rotation axis is the same as the 
symmetry axis. Thus, four rotation schemes are possible; 
oblate non-collective, prolate collective, oblate collective, 
and prolate non-collective rotations, corresponding to the 
triaxiality parameter 7 = 60°, 7 = 0°, 7 = —60°, and 
7 = —120° in the Lund convention, respectively. The ax- 
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ially symmetric limit of the RPA wobbling formalism can 
be taken for the so-called non-collective rotation schemes 
with oblate or prolate deformation, namely 7 = 60° or 
7 = —120° cases. In both cases, long-lived isomers are 
observed, but the rotational bands starting from the iso- 
mers have not been observed in the oblate non-collective 
case. On the other hand, the high-i'C isomers and the 
associated rotational bands have been known for many 
years in the Hf and W region with prolate deformation. 
Making full use of the axial symmetry, the RPA formal- 
ism has been developed |2l|, |22, 123, |24( , which is capable 
of describing the rotational band based on the high-if 
state as a multi-phonon band, i.e., the precession band. 
Recently, the same kind of rotational bands built on high- 
K isomers have also been studied by means of the tilted 
axis cranking model [25|, Hy, |27], |28( ■ 

In this paper, we would like to make a link between 
the two RPA formalisms, those for the (triaxial) wob- 
bling and for the (axially symmetric) precession motions. 
Moreover, by applying the formalism to the typical nu- 
cleus 178 W, where many high-X isomers have been ob- 
served, the properties of the precession bands are stud- 
ied in detail; not only the excitation energies but also 
the B(E2) and B(M1) values. This kind of study for 
the precession band sheds a new light on understanding 
the recently observed wobbling motion. In order to ex- 
plain the limiting procedure, we review a schematic rotor 
model in Sec. 12] while in Sec. Illll thc RPA wobbling for- 
malism and the connection to the precession band in the 
axially symmetric limit are considered. The result of cal- 
culations for 178 W is presented and discussed in Sec. II VI 
Sec. is devoted to some concluding remarks. Prelimi- 
nary results for the magnetic property of the precession 
band were already reported 29] . 



II. WOBBLING AND PRECESSION IN 
SCHEMATIC ROTOR MODEL 

The macroscopic rotor model is a basic tool to study 
the nuclear collective rotation, and its high-spin proper- 
ties have been investigated within a harmonic approxima- 
tion or by including higher order effects 0, 1^3, H3 ■ 
In this section we review the consequences of the sim- 
ple rotor model according to Ref. |30|. We use h = 1 
unit throughout in this paper. The Hamiltonian of the 
simplest triaxial rotor model is given by 
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where J's are angular momentum operators in the body- 
fixed coordinate frame, and the three moments of inertia, 
Jx, Jy and J z , are generally different. We assume, for 
definitcness, the rotor describes the even-even nucleus 
(integer spins). 

Following the argument of Ref. [3(j , let us consider the 
high-spin limit, / 3> 1, and assume that the main rota- 
tion is about the x-axis, namely the yrast band is gener- 



ated by a uniform rotation about the x-axis. Then, the 
excited band at spin I can be described by the excitation 
of the wobbling phonon, 
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where a and b are the amplitude determined by the eigen- 
mode equation, [H Iot , Xl oh ] = Lu woh (I)xl oh , at each spin 
I in the harmonic approximation. The resultant eigen- 
value w wo b(/) is given by the well-known formula, 

Wwob(/) = IyJ(l/Jy - l/Jx)(l/Jz ~ 1/Jx) 

_ fT\ I {^x ~ Jy){Jx ~ Jz) ,„s 

- Wrot(/) V ' (3) 



with the rotational frequency of the main rotation, 

U X 



(4) 



It should be noted that the triaxial deformation of 
the nuclear shape is directly related to the intrinsic 
quadrupole moments, e.g. tan7 = —V2Q22/Q20, but 
does not give a definite relation between three moments 
of inertia: One has to introduce a model, e.g. the irrota- 
tional flow model, in order to relate the triaxiality param- 
eter 7 of deformation to three inertia. Actually, if the ir- 
rotational moment of inertia is assumed, J v > J x , J z for 
the positive 7 shape, and then the wobbling frequency 
becomes imaginary. In our microscopic RPA calculation, 
the quasiparticle alignments contribute to the J x inertia 
and, as a result, the wobbling mode appears as a real 
mode even if the mean-field deformation has positive 7, 
see Refs. [H03 for details. 

The spectra of the rotor near the yrast line are given, 
in the harmonic approximation, by 

Elot {I ' n) = I ~^2j7 L + Wwob (/) ( n + 1) ' (5) 

and are composed of two sequences, the AI = 2 horizon- 
tal ones, 



£i hor) (I) =E Iot (I,n), I = n,n + 2,n + A, 



(6) 



with given phonon numbers n = 0, 1, 2, and the AI = 
1 vertical ones, 

E% et \l) = E Iot (I,I-I ), 1 = 70,10 + 1,10 + 2,... (7) 

with given band head spins lo = 0, 2, 4, both of which 
are connected by E2 transitions. The horizontal ones 
are conventional rotational bands with transition energies 
E 7 w 2ui Iot , and the AI — —2 in-band B(E2) values are 
proportional to the square of the quadrupole moment 
about the x-axis. The vertical ones look like phonon 
bands with transition energies E y w (w W ob + w ro t) and 
the AI = —1 vertical B(E2) values are 0(1/1) smaller 
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than the horizontal ones. These features are summarized 
schematically in Fig.^ In fact, the AI = —1 out-of-band 
transition was crucial to identify the wobbling motion in 
Lu isotopes [lj. If the wobbling-phonon energy Lu wo \y(I) 
is larger than the AI — 2 rotational energy AE Iot (I) — 
E xo t{I+l,n)-E Tot (I-l,n) = (21+1) /J x , both the AI = 
±1 transitions are possible. The AI = — 1 transition is 
much stronger than the AI = +1 one for the positive 
7 shape, and vice versa for the negative 7 shape, which 
also supports that the TSD bands in the Lu region have 
positive 7 shape. 
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FIG. 1: Schematic figure depicting the rotational spectra of 
a triaxial rotor Hamiltonian. The horizontal rotational bands 
are connected by solid lines, while the vertical phonon bands 
by dotted lines. 



spin of the starting state Iq = K is composed of single- 
particle alignments, i.e., the band head state is a high- if 
isomeric state. The precession band is this remaining ver- 
tical phonon band, which is actually based on a collective 
rotation about the perpendicular axis, with large angular 
momentum K along the symmetry axis. The rotational 
energy of the band is 

E high _ K (I) = -±-[I(I+l)-K% (9) 

which can be rewritten, by putting I = K + n, as 

1 n(n + 1) s 
+ 2 + 



K 



(10) 



with 



K 
J±' 



(11) 



leading to a harmonic phonon band structure with a 
one-phonon energy w prcc , when K is sufficiently large. 
This precession phonon energy coincides with the ver- 
tical transition energy given in Eq. © at the band 
head I = K. The spectra in this limit is drawn in 
Fig. El The harmonic picture holds not only for the 
energy spectra but also for the B(E2) values; for ex- 
ample, by using B(E2) cx {IfK20\IiK) 2 , one finds, in 
the leading order, B(E2;n — > n — 1) oc 3(n/K) and 
B(E2;n^ n-2) cx (3/2)(n(n-l)/K 2 ), where n = I-K 
is the number of the precession phonon quanta, so that 
the two-phonon transition is prohibited when K is large. 



Now, let us consider the axially symmetric limit |22| . 
From symmetry argument, without specifying the 7- 
dependence of three inertia, if the system is axially sym- 
metric about one of the principal axes, then two inertia 
about other two axes should coincide. Since the main 
rotation takes place about the x-axis, there are two qual- 
itatively different cases; axially symmetric about the x- 
axis or about the other y, z-axes. In the latter case, the 
main rotation is about the axis perpendicular to the sym- 
metry axis like in the case of the ground state rotational 
band, and then, apparently, the wobbling frequency Q 
vanishes; i.e., no wobbling motion takes place. In con- 
trast, if the system is axially symmetric about the main 
rotation axis x, J y = J z = J±, then the wobbling fre- 
quency © reads 



E 



^wob(I) = TP ^rotCO- 

J_L 



(8) 



Since the vertical transition energy is w w W ob + ^rot > this 
result means that the slope of the vertical sequence is 
given by I /J±, while the slope of the horizontal one by 
w ro t = 1 1 Jx- 

In reality, because the rotor model describes the col- 
lective rotation restoring the broken symmetry, J x — * 
(no collective rotation) if the x-axis is a symmetry axis. 
Thus, each horizontal band shrinks to one state, leaving 
one vertical band, whose transition energy is I/J±. The 



Precession 




® high-iT band head states 



FIG. 2: Schematic figure depicting the precession bands 
excited on high-if isomeric states. Note that the whole 
AI — ±2 horizontal sequences shown in Fig. shrink to one 
such band in the axially symmetric limit. 



III. AXIALLY SYMMETRIC LIMIT OF RPA 
WOBBLING FORMALISM 

A. Minimal coupling and RPA wobbling equation 

Microscopic RPA theories for nuclear wobbling motion 
have been developed in Refs. [Til ITU lri| . The most im- 
portant among them is that of Marshalek , where the 
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transformation to the principal axis frame (body-fixed 
frame) is performed and the theory is formulated in that 
frame. Moreover, it is shown that the RPA equation for 
the wobbling mode can be cast into the same form as 
Eq. Q if three moments of inertia are replaced with 
those appropriately defined in the microscopic frame- 
work; we call this equation the wobbling form equation. 
The adopted microscopic Hamiltonian in Ref. |14| is com- 
posed of the spherical mean-field and the quadrupole- 
quadrupole interaction (with the monopole pairing if nec- 
essary). In Ref. however, it was pointed out that the 
RPA equation could not be reduced to the wobbling form 
equation, if a most general residual interaction is used. A 
closer look into the argument in Ref. 0] shows, however, 
that the following "minimal coupling" , being used as a 
residual interaction, leads to the wobbling form equation 
as the RPA dispersion equation. 

In Marshalek's theory the rotational Nambu-Goldstonc 
(NG) modes (or spurious modes as conventionally called), 
J y and J z , play a crucial role. The RPA guarantees 
the decoupling of these modes if the selfconsistency of 
the mean-field is satisfied in the Hartree-Fock sense. In 
many cases, however, non-selfconsistent mean-fields are 
necessary; for example, the deformation is more prop- 
erly determined by the Strutinsky procedure than by 
the Hartree-Fock calculations with simple interactions, or 
one wants to study the system by hypothetically chang- 
ing the mean-field parameters, as has been done in our 
previous calculations for the nuclear wobbling mo- 
tions. Thus, we consider that the mean-field, h, rather 
than the interaction is given, and look for the residual 
interaction, Htw , which fulfills the decoupling condition 
of the NG modes within the RPA [34J . The same idea has 
been formulated in the context of the particle-vibration 
coupling theory [30l | , where the rotational invariance is re- 
stored by considering the coupling resulting from a small 
rotation about either the x, y or z-axis. Thus the mini- 
mum requirement is what we call the "minimal coupling" 
given by 



k,l—x,y,z 



XklFkFl- 



(12) 



Here, the Hermitian operator Fk and the 3x3 symmetric 
force-strength matrix \ki ar e defined as 



F k = i[h, J fe ], 

Or 1 )*! = -($\[[h,Jk],Ji]\$). 



(13) 
(14) 



with the mean-field vacuum state |4>) (Slater determinant 
if no pairing is included), on which RPA eigen-modes 
are created. If the mean-field is given by the unisotropic 
harmonic oscillator potential, the minimal coupling leads 
to the doubly-stretched Q"Q" interaction combined with 
the Landau prescription Illlailailiil^lalllSliS- 
One has to include the monopole pairing interaction in 
realistic calculations. It should be stressed that the min- 
imal coupling can be used for any type of mean-fields, 
e.g. the Woods-Saxon potential. 



For the wobbling modes in the yrast region, the mean- 
field vacuum state |#(w rot )) is obtained as the lowest 
eigen-state of the cranked mean-field Hamiltonian, 



h — h — Lo m tJ x 



(15) 



as a function of the rotational frequency uj ro t- Assum- 
ing the signature symmetry (with respect to a 7r-rotation 
about the a;-axis) of the mean-field and the conventional 
phase convention that the matrix elements of the single- 
particle operators iJ y and J z are real in the mean-field 
basis, it can be shown that the force-strength matrix \ki 
is diagonal. The excitation of the wobbling phonon cor- 
responds to the vertical AI = ±1 transitions in Sec. ITTI 
therefore only the part of the RPA equations which trans- 
fer the signature quantum number by a = 1 is relevant; 

i.e., only k,l = y,z parts of ifres in Eq. I|12fl contribute. 
It is now straightforward to follow the same procedure as 
has been done in Ref. , but with modification that the 
quadrupole field of the interaction is replaced with in 
Eq. I|12|) : Then one finds that the same RPA dispersion 
equation can be derived: 



(- 2 " ^ ot ) 



Ay(w) B z (lo) 
B y (u) A z (u) 



0, 



(16) 



where 



Ay{uj) = I - LJrotJy(uj) + UjJ yz (uj), 
A z (u) = I - LO rot J z (uj) + UjJ yz (uj), 
B y (uo) = L0j y (uj) - UJrotJyz(v), 
B z (uj) = L0J Z (UJ) ~ LL> Io tJyz(0j), 



(17) 



with the following definitions; 



I = (<P(uj Iot )\J x \<P(uj m t)) 

= 2 Jy(v)Jz{v), 

^ 2E^J y {nv) 2 

fi<u ^ v v ; 

^ 2E^J Z (^) 2 
_ ^ 2ujJ y {iiv)J z {iiv) 



(18) 



In these expressions, oj is the phonon excitation energy, 
while Epv = E^ + E v are two-quasiparticle energies with 
a = I, and J y {p.v) = (iiv\iJ y \$) (J z {p.v) = {^v\Jz\^)) 
are two-quasiparticle matrix elements of the operator 
ij y (J z ), which are associated with the vacuum state 
l^(wrot)) an d determined by the mean-field Hamiltonian 
h! in the rotating frame. It is now clear that once the 
mean-field Hamiltonian is given and the vacuum state 
l^(wrot)) is obtained, the RPA eigen-modes can be cal- 
culated without any ambiguity: This is precisely the con- 
sequence of the minimal coupling given by Eq. 1)12(1. 
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The rotational NG mode appears as a decoupled w 
uj m t solution in the RPA dispersion equation l|ltj|) ; 



rt 



1 



(iJ s + J z ) 



1 



RPA 



(tJ_) 



RPA, 



(19) 



J± — J y ± iJ z , (x-axis quantization), (20) 

where the subscript rpa means the two-quasiparticle 
transfer part (the particle-hole part if no pairing is in- 
cluded) of the operator. Note that it is normalizablc, 
[r,rt] RPA = 1, because ($\[J z ,iJ y }\$) = ($\J X \$) =1 + 
0. The cranked mean- field l|15|l describes the rotating 
state, which has an angular momentum vector aligned 
with the x-axis, and this NG mode plays a role to tilt 
the whole system by changing the x component of the 
angular momentum by —1 unit. The reason why the 
NG mode has a finite excitation energy is that there is a 
cranking term in the hamiltonian l|15l) (the Higgs mech- 
anism) . 

Finally, it has been shown by Marshalek that the 
non-NG part of the RPA dispersion equation (|16|l is re- 
duced to the wobbling form, 



J y (eff) Mji eff) M 



(iof = (c rot ) 2 



ji eff) M) 



(21) 



if three moments of inertia are replaced with microscop- 
ically defined ones in the following way; 



x 



i = (^KoOI^Kot)) 

^rot W ro t 



jj eff >(u;) = J„(u) - J vz {oj) 



B z {u)' 
A z (co) 



(22) 



Since the y- and z-effective inertia are w-dependent, the 
equation is non-linear and they are determined only after 
solving it. 

As for the electromagnetic transition probabilities, 
Marshalek proposed a 1/ /-expansion technique by utiliz- 
ing the perturbative boson expansion method hj. The 
Al = =pl E2 and Ml vertical transitions from the one- 
phonon wobbling band to the yrast band, discussed in 
Sec. m can be calculated within the RPA, which is the 
lowest order in 1/J, as 

B(E2;I±1 - I) « M[Q2 T uXl oh }\<P)\ 2 , (23) 
B{M1; I ± 1 - I) » | (4>\ [ MlTl , Xl h ]\<P)\ 2 , (24) 

where A^ ob is the wobbling phonon creation operator, 
and the E2 and Ml operators quantized with respect to 
the x-axis, 



Q2±l - ^T^Qzi ± Qi2)i 

Ml±l = ± ~^( i ^y T 



(25) 
(26) 



are introduced (see also Ref. H^H^)- Here (K = 

0,1,2) are electric quadrupole operators (z-axis quanti- 
zation) with a good signature, 



?22 




a=l 
Z 



E(rf } > 



(27) 



and /ifc (k = x, y, z) are magnetic dipole operators, 
li k = \ — /ijv J2(9 ( i T) h + 9^ T) s k ) a , (r = tt, v). (28) 



B. Axially symmetric limit and RPA precession 
equation 

If the deformation is axially symmetric about the x- 
axis, the angular momentum is generated not by the col- 
lective rotation, but by the alignment of the angular mo- 
menta of quasiparticles along the symmetry axis. The 
mean-field vacuum state |#), a high- if state, is a multi- 
ple quasiparticle excited state, and its spin value is the 
sum of the projections, f2^, of their angular momenta on 

the symmetry axis; I = K = X^° CC ^ ^p-i i- e -, the time re- 
versal invariance is spontaneously broken in |<P). In this 
case, the cranking term in Eq. (|15l) does not change the 
vacuum state |#), so that the rotational frequency uj mt is 
a redundant variable: All observables should not depend 
on Lu mt . It is reflected in the fact that the quasiparticle 
energies linearly depend on the rotational frequency: 



Efj,(u} TO t) = E° — W ro t^p 



(29) 



where E® are quasiparticle energies for the non-cranked 
mean- field Hamiltonian h. Since the cigcn- value of J x , 
f2, is a good quantum number, it is convenient to rewrite 
the RPA dispersion equation l|l(j|) in terms of the matrix 
elements of J± rather than iJ y and J z : After a little 
algebra, the equation decouples into two equations, 



(u ± CJ rot ) S±i(u> ± CJ r ot) = 0, 



(30) 



where the functions S p (uj) with p = ±1 determine the 
Afl = ±1 solutions, respectively, and are given by 

1 ^ f (E^y ± u rot )\J±{nv)\ 2 



\1<V 



E^v ± ijj ro t — u) 
(Ejw =F ^rot)| J t (a^)| ; 

En,v T ^rot + W 



(31) 



The precession is a AQ = +1 mode, as is clear from 
the rotor model in Sec. [HJ and then only the Al = — 1 
E2 and Ml transitions are allowed; i.e., their Al = +1 
probabilities vanish in Eqs. (|23|l and l|24|) because the 
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two RPA transition amplitudes, (#| [Q 21 \ -^wob] 1^) an< ^ 

(^HQ^^^wob]!^)' are the same m their absolute value 
with the opposite sign; a corresponding relation holds for 
the Ml amplitudes. 

On the other hand, the y- and z-inertia are the same 
due to the axial-symmetry about the s-axis, and then, 
just like Eq. JBJl, Eq. f2*TJ) reduces to 



dispersion equation and vice versa. For the physical 
AQ = +1 modes, the eigen-energies of the wobbling dis- 
persion equation (plOf and the precession one lf3f)|l are 
related as 



<^wob — k^prec ^rot • 



By comparing it with Eq. H32[l . we obtain 



LU = ± 



K 



Tw Iot (Af2 = ±l) 



(32) 



K 



^prcc 



off) ' 



where I=($\J X \$) is denoted by K, and the perpendic- w ith J r j cff) being written as 



ular inertia J^± S \u)) 
written as 



J y (eff V) = Jr\u) is simply 



(cff) 



j[ eff) H = Jj_(u)tJ V z{u) (AQ = ±1), (33) 

With J ± (U) = Jy(u) = J Z {UJ). 

The vibrational treatment of the rotational band built 
on the high-if isomeric state in terms of the RPA has 
been done for a harmonic oscillator model in Refs. 0, 
122^ . and for realistic nuclei by employing the Nilsson 
potential in Ref. |2^, followed by calculations with the 
Woods-Saxon potential in Ref. [24| . The residual interac- 
tion adopted in Refs. 13 is derived by ap plyi ng the 
vibrating potential model of Bohr-Mottelson |30j to an 
infinitesimal rotation about the perpendicular axis, and 
equivalent to the minimal coupling (|12J) : In the axially 
symmetric case, 



H i3at = ~K(FlF + + FlF^), 
with F± being defined by using J± in Eq. 



(34) 



F± = -[h,J±], k 

K 



■ : {$\[[h,J-],J+]\*). (35) 



Note that the mean- field state \<P) is now a multi- 
quasiparticle excited state for the non-cranked mean-field 
Hamiltonian h, and so w ro t does not appear, although 
it can be used as the "sloping Fermi surface" to obtain 
optimal states 0]: The cranking procedure is totally 
unnecessary in this approach. 

The resultant RPA dispersion equations are given for 
the parts associated with the fields F± separately; 

wS±i(w)=0, (36) 

where the functions S±i(u) are defined by 



1 J2 \ e%\j ± (H\ 2 



E° 



(37) 

which turn out to be the same functions as Eq. I|31|l be- 
cause of the property (|29() of quasiparticle energies in 
the non-collective rotation scheme. It is worth mention- 
ing S + i(uj) = —S^i(—lj), so that Ail = —1 modes are 
obtained as negative energy solutions of the AQ = +1 



J 



(off) 



5D 



\J + {H\' 



|j_(/«/)| s 



(38) 



(39) 



(40) 



which is the microscopic RPA version of Eq. in 
Sec. [n] This ,7^ does not depend on w rot , while both 
J± = J y = J z and J yz in Eq. do. This result can 
also be obtained directly from the precession dispersion 
Eq. H36[) . Note that the perpendicular inertia (|40|l re- 
duces to the Inglis cranking inertia (or that of Belyaev if 



pairing is included) in the adiabatic limit w p 



0. 



The reason why the a> ro t-dependent wobbling eigen- 
energy and the w ro t-independent precession eigen-energy 
is related in a simple way l|38|l is that the RPA treatment 
in Refs. pll I22I I23L l24| is formulated in the laboratory 
frame, while Marshalek's wobbling theory in the princi- 
pal axis frame (body- fixed frame). The energies in the 
laboratory frame, E^ , and in the uniformly rotating 
frame described by the cranked mean-field, E^ VR \ are 
related by £( UR ) = E^ ~ Qu) ro t for the state which has 
a projection, Q, of angular momentum on the rotation 
axis. Moreover, the energies in the principal axis and the 
uniformly rotating frames are the same under the small 
amplitude approximation in the RPA. Thus the differ- 
ence of phonon energies in l|38|l comes from the difference 
of coordinate frames where the two approaches are for- 
mulated. The rotational NG mode l|19|) . appears at 
zero energy in the precession dispersion Eq. (|36[) by the 
same reason. The transformation between the laboratory 
and the principal axis frames have been discussed more 
thoroughly in Refs. fli] and ^2}. 

As for the electromagnetic transition probabilities in 
the precession formalism |2J, [24j , the RPA vacuum state 
I RPA) is considered to be a stretched eigen-state of 
the angular momentum, \I = K,M = K), because 
r|RPA) = for the NG mode (F oc (J + ) RPA ). In 
the same way, the AQ = +1 one-phonon precession state 
A^JRPA) corresponds to \I = K + 1, M = K + 1), be- 
cause rAt roc |RPA) = [r,At roc ]|RPA) = 0. Then, by 
using the Wigner-Eckart theorem, we obtain, for exam- 
ple, 



(I = K\\M{E2)\\I = K + 1) = V2K+1 
( RPA I Q 2 - 1 A^ rec I RPA) 
X (K+ IK+ 12 -1\KK)' 



(41) 
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Thus, by inserting explicit expressions of the Clebsch- 
Gordan coefficients, one finds 



B(E2;K + 1 -> K) 



K + 2 
K 



|<*|[Q a _i,JTt p)f 



(42) 

S(Ml;Jf + 1 -> X) = K^^X^Jl^^ (43) 

which coincide, within the lowest order in 1/K, with 
Eqs. (|23H and l|24|) in the wobbling formalism. 



IV. RESULT AND DISCUSSION 

A. Calculation of precession bands in 178 W 

In the previous papers 0, we have studied the 
wobbling motions in the triaxial superdeformed bands in 
Hf and Lu isotopes. As it is demonstrated in the previ- 
ous section, the precession mode can be described as an 
axially symmetric limit of the wobbling formalism. Thus 
we have performed calculations of the precession bands 
in 178 W, for which richest experimental information is 



available 43, 44, 45]. Exactly the same wobbling formal- 
ism is used, but taking the prolate non-collective limit 
suitable for high- if isomers, i.e., the triaxiality parame- 
ter 7 = —120° in the Lund convention. The first result 
for this nucleus, concentratin g on the magnetic property, 
was reported already in Ref. |29( . 

The procedure of the calculation is the same as in 
Refs. [ll El E|: The standard Nilsson potential 
is employed as a mean-field with the monopole pairing 
being included; 

h = h ms {e 2 ,f)- At(P} + Pt)- J2 X - N - ( 44 ) 



Here the €4 deformation is neglected and all the mean- 
field parameters are fixed for simplicity. There are a few 
refinements of calculation, however: 1) the difference of 
the oscillator frequencies for neutrons and protons in the 
Nilsson potential is taken into account, and the correct 
electric quadrupole operator is used, while Z/A times 
the mass quadrupole operator was used previously, 2) the 
model space is fully enlarged; N osc =3— 8 for neutrons and 
2—7 for protons, which guarantees the NG mode decou- 
pling with sufficient accuracy in numerical calculations. 
As for the point 1), usually, QW « (Z/A){Q^ + Q«) 
holds for static and RPA transitional quadrupole mo- 
ments in stable nuclei, and therefore, the simplification 
in the previous paper was a good approximation. It is, 
however, found that QW is appreciably smaller, by about 
4-8%, than (Z/A)(Q^ + Q«) in 178 W. Thus, in this 
paper, we make a more precise calculation using the elec- 
tric (proton) part of the quadrupole operator. 

The calculation is performed for the high-if isomeric 
configurations listed in Table [I] they cover almost all 
the multi-quasiparticle states higher than or equal to 
four (more than or equal to two quasineutrons and two 



quasiprotons) , on which rotational bands are observed. 
The quadrupole deformation is chosen to be e 2 = 0.240, 
which reproduces in a rough average the value Qq = 7.0 
b for the configurations in Table [I] assumed in the experi- 
mental analyses |44l l45j . The pairing gap parameters are 
taken, for simplicity, to be 0.5 MeV for two-quasiparticle 
configurations, and 0.01 MeV for those with more than 
or equal to four-quasiparticles, both for neutrons and 
protons. Chemical potentials A r (r = v, it) are always 
adjusted so as to give correct neutron and proton num- 
bers. These mean that the choice of parameters in this 
work is semi-quantitative. As is explained in detail in 
Sec. IIIII the final results do not depend on the cranking 
frequency w rot at all for the non-collective rotation about 
the x-axis. We have confirmed this fact numerically and 
used w ro t = 0.001 MeV in actual calculations (Note that 
the wobbling RPA formalism requires a finite rotational 
frequency in numerical calculations) . No effective charge 

is used for the E2 transitions, and gi eS ^ — 0.7gi hcc ^ is 
used for the Ml transitions as usual. 
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FIG. 3: Dependence of numerical results on the pairing gap 
parameter A = A n = A v . In the upper panel, the excitation 
energies are shown, while in the lower panel are shown the 
RPA transition amplitudes for the electric Q22 operator 1271 . 
The solid curves show the results for the precession modes 
excited on the 25 + and 30 + high-Tf states in 178 W, and the 
dotted, dashed and dot-dashed curves represent those for the 
7-vibrations in 166 Er, 168 Yb and 178 Hf, respectively. 



We have checked the dependences of the results on the 
variations of the deformation parameter e 2 and pairing 
gaps. Those on the pairing gaps are shown in Fig. [3] 
In this figure, the excitation energy U) and the RPA 
transition amplitude for the electric Q 22 ' operator l|27|l , 
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TABLE I: Configurations assigned for high-Jf isomers in 178 W |43l I44I |45| . which are used in the RPA calculations for 
the precession bands excited on them. The experimental values of the precession one-phonon energy, ui^i c = Ek(I = K + 
1) - Ek{I = K), are also tabulated in the last column. The neutron states are 1/2" [521], 5/2" [512], 7/2" [514], 7/2+ [633], 
9/2+ [624], and 7/2" a [503]. The proton states are 1/2"[541], 5/2+[402], 7/2+[404], 9/2~[514], and ll/2~[505]. The bold 
letters indicate the h 9 / 2 proton and the 1-12./2 neutron quasiparticles. 



K 




Neutron configuration 




Proton configuration 


, ,cxp f\ rr 
^prcc 


13" 






7 /*>+ 7 ! r )~ 




5/2 > 7/2^ 


164 


1 A + 

14 






7/2+, 7/2" 




5/2 , 9/2 


150 


15+ 






7/2+, 7/2" 




7/2+, 9/2- 


207 


18~ 






7/2+, 7/2- 




1/2-, 5/2+, 7/2+, 9/2- 


184 


21" 




5/2" 


, 7/2+, 7/2", 9/2+ 




5/2+, 9/2- 


362 


22" 




5/2" 


, 7/2+, 7/2", 9/2+ 




7/2+, 9/2- 


373 


25+ 




5/2 _ 


, 7/2+, 7/2-, 9/2+ 




1/2-, 5/2+, 7/2+, 9/2- 


288 


28" 




5/2 _ 


, 7/2+, 7/2", 9/2+ 




1/2", 7/2+, 9/2", 11/2- 


328 


29+ 


5/2" 


, 7/2+ 


, 7/2-, 9/2+, 1/2- . 


, 7/2-" 


1/2-, 5/2+, 7/2+, 9/2- 


437 


30+ 




5/2- 


, 7/2+, 7/2-, 9/2+ 




5/2+, 7/2+, 9/2", 11/2- 


559 


34+ 


5/2" 


, 7/2+ 


, 7/2-, 9/2+, 1/2- . 


, 7/2" a 


5/2+, 7/2+, 9/2", 11/2" 


621 



TABLE II: Mean-field parameters used in the calculation 
for the 7- vibrations on the ground states (7 = 0°), and ob- 
served excitation energies of 7- vibrations |47| . The £2 values 
are taken from Ref. |4q| . where they are deduced from the 
measured B(E2 : 0+ — > 2+) values. The even-odd mass dif- 
ferences are calculated by the third order difference formula 
with using the binding energy data in Ref. |49|. 



Nucleus 


£2 


Z^(MeV) 


zW(MeV) 


o4 xp (MeV) 


166 Er 


0.272 


0.966 


0.877 


0.786 


168 Yb 


0.258 


1.039 


0.983 


0.984 


178 Hf 


0.227 


0.694 


0.824 


1.175 



Q = \([Q 2 2 1 -^prcc]) 1 1 which is a measure of the E2 collec- 
tivity, for the precession modes excited on the K = 25+ 
and K = 30+ configurations, are shown as functions 
of the pairing gap, A = A„ = A v (the common value 
for protons and neutrons). For reference sake, are also 
included the results for the 7-vibrations on the ground 
states, i.e., the AK = ±2 vibrational mode excited on 
the 7 = 0° prolate mean- field (without cranking), for 
166 Er, 168 Yb and 178 Hf nuclei. Note that the meaning of 
the operator is different for 7 = —120° and 7 = 0° 
shapes, so that the comparison of the magnitude of the 
amplitude Q is not meaningful between the precession 
mode and the 7-vibrational mode. As it is stressed in 
Sec. IIII Al the precession mode is calculated without any 
ambiguity once the mean- field is fixed; we have just used 
the same parameters explained above with an exception 
that the pairing gaps are varied. The situation for the 
7-vibration is different; one has to include components 
other than the minimal coupling, l|12l) or l|34l) . We have 
used the K = 2 part of the doubly-stretched Q"Q" force, 
and the force strength is determined in such a way that 
the calculations with adopting the even-odd mass differ- 



ences as pairing gap parameters reproduce the experi- 
mental energies of the 7-vibration; see Table [D] for the 
parameters and data used. Then, with the use of the 
force strength thus fixed, calculations are performed with 
varying the pairing gaps. 

As is clearly seen in Fig. [21 the reduction of pair- 
ing gaps makes the excitation energies of 7-vibration 
change in various ways depending on the shell struc- 
ture near the Fermi surface; i.e., the distribution of 
the Afl = ±2 quasiparticle excitations, which have 
large quadrupole matrix elements. The energy becomes 
smaller and smaller in the case of 166 Er, and finally 
leads to an instability (w 7 — > 0); accordingly the tran- 
sition amplitude Q diverges. No instability takes place 
in the case of 168 Yb and the excitation energy decreases 
with decreasing A, while it is almost constant for the 7- 
vibration in 178 Hf. However, the transition amplitudes Q 
reduce by about 40—60% with decreasing A except for 
166 Er. These are well-known features for the low-lying 
collective vibrations; namely, the collectivities of the vi- 
brational mode are sensitive to the pairing correlations; 
especially enhanced by them. In contrast, for the case 
of the precession modes, the excitation energies are sta- 
ble and transition amplitudes are surprisingly constant 
against the change of the pairing gap. This is a feature 
common to the wobbling mode excited on the triaxial su- 
perdeformed band 01 • Although both the precession (or 
the wobbling) and the 7-vibration are treated as vibra- 
tional modes in the RPA, the structures of their vacua 
are quite different; the time reversal invariance is kept 
in the ground state while it is spontaneously broken in 
the high-spin intrinsic states. Since the precession or 
the wobbling is a part of rotational degrees of freedom, 
this symmetry-breaking may be an important factor to 
generate these modes. It should be mentioned that the 
transition amplitude Q for 166 Er leads to about a factor 
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of two larger B(E2 : 2+ — > 0+) value than the observed 
one in the present calculation, in which the model space 
employed is large enough. The RPA calculation overesti- 
mates the B(E2) transition probability for the low-lying 
7- vibration if are used the Nilsson potential as a mean- 
field and the simple pairing plus Q"Q" force as a residual 
interaction [50| . 

There are many RPA solutions in general, and it is 
not always guaranteed that the collective solution ex- 
ists. In some cases collective solutions split into two or 
more, whose energies are close, and the collectivity is 
fragmented (the Landau damping), or the character of 
the collective solution is exchanged. Moreover, in the 
case of precession-like solutions, the fl = ±1 modes in- 
teract with each other, as it was shown in Ref. |23j |. In 
fact, when the deformation is changed, it is found that 
the precession mode on the K 11 = 15 + configuration dis- 
appears for £2 > 0.250, and that on the = 14~ splits 
into two for e 2 > 0.245. Similar situations also occur 
when changing the pairing gap parameters in a few cases. 
Apart from these changes, the results are rather stable 
against the change of the mean-field parameters. The 
fact that we have been able to obtain collective solutions 
for all the cases listed in Table [I] indicates that our choice 
of mean-field parameters are reasonable if not the best. 

hu>(MeV) 
0.7 




state is given. With a closer look, however, one finds 
deviations especially at K* = 18~, 25 + , 28~, and 29 + : 
The precession energies on them are smaller in compar- 
ison with others, but the calculated ones are too small. 
Low calculated energies correspond to large perpendicu- 
lar moments of inertia as is clearly seen in Fig. 03 These 
four configurations contain the proton high-j decoupled 
orbital (i.e. with Q = ±1/2) 7r[541]l/2 _ originating from 
the hg/2, whose decoupling parameter is large. Occupa- 
tion of such an orbital makes the Inglis moment of in- 
ertia, which is given by Eq. I|4l)|l with setting w prcc = 0, 
diverge due to the zero-energy excitation from an occu- 
pied i? = +1/2 quasiparticle state to an empty —1/2 
state. The reason of too large moment of inertia may be 
overestimation of this effect for the proton contribution 
in the calculation. The large effect of this Tthg/2 orbital 
on the moment of inertia has been pointed out also in 
Refs. El EEL 



J{h 2 /MeV) 




20 25 30 35 

K(h) 



FIG. 5: The moment of inertia associated with the precession 
bands built on high-if configurations. The RPA effective in- 
ertia il HH are shown by filled circles connected by solid lines, 
the proton part of them by filled squares connected by dot- 
ted lines, and the neutron part by filled triangles connected 
by dashed lines. The crosses are those extracted from the 
experimental spectra according to the simple relation 13911 . 



FIG. 4: Excitation energies of the one-phonon precession 
modes excited on high- K configurations. Calculated ones are 
denoted by filled circles connected by solid lines, and experi- 
mental ones by crosses. Data are taken from Refs. |44ll45||. 

Figure 0] presents the calculated and observed relative 
excitation energies of the first rotational band member, 
Ei = k+i — Ei=k, i-e., the one-phonon precession ener- 
gies. Corresponding perpendicular moments of inertia, 
Eq. (|39|l . are shown in Fig. [SJ where the contributions 
to the inertia from protons and neutrons are also dis- 
played. Our RPA calculation reproduces the observed 
trend rather well in a wide range of isomeric configura- 
tions, from four- to ten-quasiparticle excitations. This is 
highly non-trivial because, as was stressed in Sec. lIIII we 
have no adjustable parameter in the RPA for the calcula- 
tion of the precession modes once the mean-field vacuum 



Except for the case of four configurations including 
the 7r[541]l/2~ orbital, the values of moment of iner- 
tia are about 50—80 fi, 2 /MeV, which are smaller than 
the rigid-body value, J rig =87.8 h 2 /MeV, and consid- 
erably larger than the ground state value, Jg r =28.3 
fi 2 /MeV. Here J r - lg is calculated by assuming the 178 W 
nucleus as an ellipsoidal body with e 2 — 0.240 and 
ro = 1.2 fm, and J ET by 3/E 2 +. The pairing gaps are 
already quenched in the calculation for more than or 
equal to eight-quasiparticles (four-quasiprotons and four- 
quasineutrons) configurations (K > 25 + ). The value 0.5 
MeV of the pairing gap used for two-quasiparticles con- 
figurations is already small enough to make the moment 
of inertia quite large. It is also noticed that the mo- 
ment of inertia decreases with increasing K , which is op- 
posite to the intuition, and clearly indicates the impor- 
tance of the shell effect for the moment of inertia 1531. 
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In Refs. |25l l44j. the angular momentum of the preces- 
sion band is divided into the collective and aligned ones; 
the inertia defined in Eq. 1)39(1 includes both of them. It 
is shown that the collective inertia, in which the effect 
of the aligned angular momentum of the high-j decou- 
pled orbital is removed, takes the value 50—60 ft 2 /MeV 
consistent with the other configurations. As is shown 
in Fig. [S] the proton contribution to the inertia is about 
20—30% (except for the four configurations above), which 
is considerably smaller than Z/A, but consistent with the 
calculated value for the <?^-factor in the ground state ro- 
tational band (see below). 

As for the electromagnetic transitions in the rotational 
bands built on high- A" isomers, the strong coupling rota- 
tional model |3pJ is utilized as a good description. The 
expressions for B(E2) and B(M1) are well-known: 



B(E2 : I = A + 1 -> K) TO 



16tt 



e 2 Q 2 (A + 1 A20|AA) ; 



15 1 
16tt K 



e Qq, 



(45) 
(46) 



B(Ml : I = K + 1 -> A) rot 

= (9K ~ 9R? K 2 (A + 1 A 10| AA) 2 (47) 



-^/4f (9k - 9nf A, 



(48) 



where, in the last lines, the Clebsch-Gordan coefficients 
are replaced with their lowest order expressions in 1/A'. 
Qq and (gx — go) can be extracted from experiments; 
the sign of the mixing ratio is necessary to determine 
the relative sign of them. These quantities are calculated 
within the mean-field approximation, 



Qo = V~(Q 20 ) = (£(2* 



/47T (fjt x ) 

9K= \l^r 



2 y 2 -z 2 t ] 



3 hn{Jx 



(49) 

9r = \-z rn—> ( - 50 ) 

V 3 HN\Jx)gi 



where ( ) means that the expectation value is taken 
with respect to the high- A configuration (7 = —120°), 
e.g. (J x ) = A, while ( ) gr with respect to the ground 
state rotational band (7 = 0°). The latter expectation 
value is calculated by the cranking prescription l|15|l . with 
the same 62, and with the even-odd mass differences as 
pairing gaps. The value of gn is thus w rot -dependent, 
but its dependence is weak at low frequencies, so that we 
take the value gn = 0.227 obtained at LU TOt — > 0, which is 
much smaller than the standard value, Z/A — 0.416. 

On the other hand, B(E2) and B(M1) are calculated 
by Eqs. i|23|l and (I24II . respectively, in the wobbling for- 
malism which is in the lowest order in 1/A. By equating 
these expressions with those of the rotational model, 146|) 



and H48|) . we define the corresponding quantities in the 
RPA formalism by (A = (J x )) 



(q ; 



RPA 



'WnK 1 

15 e 



Xl iec ,Q 2 _ 1 ]), (51) 



(9K - 9r)kPA = \/xT7 ([^prec)Ml-l]}- (52) 

Only their relative phase is meaningful, and the over- 
all phase is chosen in such a way that (Qo)rpa is pos- 
itive. We compare calculated values of Qo m the usual 
mean-field approximation l|49|) and in the RPA formal- 
ism (|51|l in Fig. for all high-A configurations listed 
in Table [I] These two calculated Qo's roughly coin- 
cide with each other, but appreciable deviations are seen 
for the A T = 18", 25+, 28", and 29+ isomers: The 
high-j decoupled orbital 7r[541]l/2" has a large prolate 
quadrupole moment, so that its occupation generally 
leads to a larger value of Qq. This is clearly seen in 
Fig. 0even if e 2 is fixed in our calculation. See Ref. [53| . 
for example, for the polarization effect of this high-j or- 
bital on Qq. It is, however, noticed that the effect is even 
larger in the RPA calculation just like in the case of the 
excitation energy in Fig. For the 34 + isomer we have 
found a less collective RPA solution at a lower energy, 560 
keV, which has about 80% of the (Qo)rpa value of the 
most collective one presented in the figure. The reason 
why (Qo)rpa for the 34 + isomer is considerably small is 
traced back to this fragmentation of the precession mode 
in this particular case. This kind of fragmentation some- 
times happens in the RPA calculation. 




20 25 
K(h) 



FIG. 6: Quadrupole moments Qo for high- if configurations. 
Those calculated by the RPA, Eq. 1511 . are denoted by filled 
circles connected by solid lines, while those by the mean-field 
approximation, 149^ . by filled triangles connected by dotted 
lines. 

In Fig. we compare the effective (gx — go) factors 
extracted from the experimental data and those calcu- 
lated in two ways; Eq. ((5Tj|) and Eq. (1521) . As for the 
observed ones they were determined |44l |45j from the 
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K(h) 



suits indicate, however, that the RPA treatment of the 
rotational excitations is valid; especially it gives a reliable 
microscopic framework for studying the wobbling motion 
recently observed. 

The reason why the RPA precession mode gives the 
B(E2) and B(M\) similar to those calculated according 
to the rotational model is inferred by taking the adiabatic 
limit (wprec — > 0) of the RPA phonon creation operator. 
It has been shown in Ref. that the precession phonon 
can be explicitly written up to the first order in oj prec as 

-\J+ + ^proc l 7T 1 *6 ) +)RPA 



At 

prcc 



= (J+ + KiO+) RPA . 



(53) 



FIG. 7: Effective (gpc — <?ii)-factors for high- if configurations. 
Those calculated by the RPA, Eq. I52II . are denoted by filled 
circles connected by solid lines, while those by the mean-field 
approximation 14911 by filled triangles connected by dotted 
lines. Those extracted from the experimental data |44l |45| 
are shown by crosses. 



branching ratios of available lowest transitions in respec- 
tive rotational bands, by using the rotational model ex- 
pressions (|45fl and l|47|l with Qq = 7.0 b being assumed. 
In this way, absolute value \gx — gn\ is obtained, and we 
assume that its sign is determined by that of the calcu- 
lated E2/M1 mixing ratio in the RPA result. Accord- 
ingly, some care is necessary to compare the experimen- 
tally extracted ^-factors with calculations. The agree- 
ment between the observed and calculated ones is semi- 
quantitative, but the RPA result follows the observed 
trend rather well. Compared to the RPA ^-factors, those 
calculated by the mean-field approximation are poorer. 
Again, the two calculations deviate appreciably for the 
K* = 18", 25+, 28", and 29+ configurations, where the 
high-j decoupled orbital ir [541] 1/2", which has a large 
positive g-factor, is occupied. The difference between 
the mean-field (j)k — 9r) and (gx — <?,r)rpa is further 
discussed in the next subsection by studying the adia- 
batic limit of the precession mode in the RPA. 



Interpretation of the result in the adiabatic 
limit 



As it is demonstrated in the previous subsection, the 
RPA calculation reproduces the precession phonon ener- 
gies without any kind of adjustments. The electromag- 
netic properties obtained through the RPA wobbling for- 
malism are in good agreement with those of the strong 
coupling rotational model, where the quadrupole mo- 
ments and the effective g-factors are calculated within 
the mean-field approximation. Since the rotational band 
is described as multi-phonon excitations in the RPA wob- 
bling (or precession) model, it is not apparent that two 
models lead to similar results for observables. Our re- 



Here the angle operator + is defined by 

0± = 0y±iO z , 



[h,i0 k 



1 



(54) 



where J" is the Inglis cranking inertia and given from 
the effective inertia (|40|) by setting uj picc = 0. These 
angle operators possess desired properties, 

([©*, Ji]) = i5 H . (55) 

For the E2 transitions, the contribution of the 0+ part 
in Eq. (|53|l is proved to be negligible, if the harmonic 
oscillator potential is taken as a mean-field; 



[^prec;Q2-l] 



= <[J+,Q 2 -l]) 



(56) 

(<5o)rpa m the adiabatic 



which precisely means Qo 
limit. 

As for the Ml transitions, however, the 0+ part also 
contributes: 



[^prec>Ml-l]) 




Mi-i 



- 77| <[Mi-i,? ; 0+]) 



which gives (g K - go) « (gK 



9r <-> 9r 



(57) 

cm) rpa if we identify 

(58) 



3 hn 



This identification is reasonable: The magnetic moment 
operator [ii-i possesses a property of angular momen- 
tum and is approximately proportional to J_ . Then the 
expectation value of the right hand side of Eq. (frTSf is 
expected to depend only weakly on the high- if configu- 
ration because of Eq. (|55|l . More precisely, if the opera- 
tors J_, 0+ and are divided into the neutron and 
proton parts like 



J- = J 



J 



0, 







(t) 







f^( 5 Wj W + 9 Mj M ), (59) 
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then the following relation is derived, 



9R 



(60) 



because of ([ji r) ,i<9i T) ]) = 2J™ [t) fjf with Jl 



Tr,v). With a cruder estimate 



([ji T \ ~ 2N T /A (r = one finds a constant 

9R ~ ^2i T N r g ^ /A, which gives a classical result, Z/A, 
by setting #W = 1 and g {v) = 0. 




FIG. 8: Comparison of deduced gn from two calculations. 
The one from qk — (qk — 5h)rpa is denoted by filled circles 



connected by solid lines, while the quantity J 



(off, it) 



/(J 



(cff.vr) 



+ 



J 



(cft>h 



by filled squares connected by dotted lines. 



An approximate relation g R = J M /(J M + jM), 
which corresponds to Eq. (|60|l with <jrM = 1 and 
g( u ' — 0, has been used for the ground state rota- 
tional band, i.e., the case of collective rotations p54j . 
It seems, however, difficult to justify a similar relation, 
g R = ZjW /(ZjW +NJW), which is used in Ref. 0. 
Thus, the "rotor g-i actor" g R is not a common constant, 
but it also depends on the high-i'T configurations as the 
intrinsic g-factor gx does. In order to see how the ap- 
proximate relation (|f>()|> holds, we compare, in Fig. |H1 
the two calculated quantities, gx — {gx — <7_r)rpa and 



J'±' V ^), where the cranking inertia 

j" yT> , which diverges when the 7r[541]l/2~ orbital is oc- 
cupied, is replaced with the neutron or proton part of 
the effective inertia l|40[) . see also Fig. [5j As is seen in 
the figure, these two quantities are in good agreement 
with each other, again, except for the K* = 18~, 25 + , 
28 _ , and 29 + configurations, where the high-j decou- 

(cff.jr) 



pled orbital is occupied and j r | cff,7r ) J{JJ 



is very large. The excitation energies are underesti- 
mated for these high-if configurations. Therefore, the 
proton moments of inertia are overestimated for them; 
in fact the proton contributions are considerably larger 
than the neutron ones in these configurations as is shown 
in Fig. [SJ Apart from these four configurations, the 



deduced (/^-factors in Fig. [S] are similar to the ground 
state value, 0.227, though it is appreciably different from 
the standard value, Z/A = 0.416. For reference, the 
cranking moment of inertia for the ground state rota- 
tional band calculated using the even-odd mass differ- 
ences as pairing gaps is = 22.7 h 2 /MeV (about 
80% of the experimental value, see the previous subsec- 
tion). The proton contribution to it is 6.1 h 2 /MeV and 
j{ cr )«/(j[ cr )W + jf r )M) = 0.269, which is slightly 
larger but consistent with the calculated ground state g R 
value, 0.227. 

The above results indicate that the rotor g R should 
be considered to depend also on the intrinsic configura- 
tions, but the dependence is conspicuous only for those 
including the high-j decoupled orbit, which has a large 
decoupling parameter as well as a large g-i actor. The 
reason why the effective (gx — <7ij)-factors of the RPA 
calculation reproduce the experimentally extracted ones 
better than those of the mean- field g-factors is inferred as 
follows. Since, as is well known, the (/-factors of proton 
orbitals are much larger than those of neutron orbitals, 
the amount of the proton contribution is overwhelming 
for the mean- value (fi x ) in comparison with that for (J x ). 
Considering this fact together with the overestimation of 
the proton moments of inertia mentioned in the previ- 
ous paragraph, it is likely that the calculated values of 
g K © for the K* = 18", 25+, 28", and 29+ configu- 
rations with a proton high-j decoupled orbital are also 
overestimated. In the mean-field calculation, the calcu- 
lated values of (gx — 9r) for those configurations are 
thus relatively large, because the common ground state 
<?,R-factor (|50J) is used. This trend can be seen also in the 
similar type of mean- field calculations in Refs. liiL lial 
(Note that different (/^-factors are used in 0] and |45|). 
In the RPA calculation, however, the rotor (/-factor is 
given by g R , 1(58(1 or 1)60(1 . which is overestimated for these 
four configurations (see Fig. [HJ - Thus, the overestimation 
of two (/-factors may largely cancel out in the resulting 
{dK ~ 9r)rpa values, yielding a reasonable agreement 
with the experimental data seen in Fig. 

The realistic mean-field is not very different from the 
harmonic oscillator potential so that the approximate 
equality for the E2 operator is expected to hold in 
general cases. However, it is not very clear to what extent 
this equality holds: It is a subtle problem whether or not 
the adiabatic approximation holds because the precession 
phonon energies are 200 to 600 keV, which are not neg- 
ligible compared to the quasiparticle excitation energies 
(Note that the pairing gap is quenched in high-iif con- 
figurations). In addition to the deviations caused by the 
non-adiabatic effects, it should also be noticed that the 
adiabatic approximation itself breaks down if one quasi- 
particle in a pair of high-j decoupling orbits (f2^ = ±1/2) 
is occupied, because the Inglis cranking moment of iner- 
tia diverges due to the zero-denominator. In such cases, 
the present RPA calculation eventually overestimates the 
moment of inertia, although it does not diverge. This ef- 
fect is also reflected in the calculated transition moments 
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(Qo)rpa and the effective g- factors, which are rather dif- 
ferent from the values given by the mean-fields. Whether 
or not the RPA calculation gives reliable results for such 
cases, where the non-adiabatic effect is large, is an im- 
portant future issue. The direct measurement of Qq (i.e. 
B(£'2)-value) for the precession band is desirable for this 
purpose. 

V. CONCLUDING REMARKS 

We have investigated the precession bands, i.e., the 
strongly-coupled rotational bands excited on high- if in- 
trinsic configurations by means of the random phase ap- 
proximation, the microscopic theory for vibrations. It is 
emphasized that this precession mode is related to the 
three dimensional motion of the angular momentum vec- 
tor in the principal axis frame (body- fixed frame), and 
can be considered to be a limiting mode of the wobbling 
motion in the triaxially deformed nucleus. It is demon- 
strated that the observed trend of the precession phonon 
energies in 178 W is well reproduced by the RPA calcula- 
tion: This is highly non-trivial because we have employed 
the minimal coupling interaction, which is determined by 
the mean-field and the vacuum state based on it, and so 
there is no adjustable force parameters whatsoever. 

The electromagnetic properties, the E2 and Ml tran- 
sition probabilities, are also important for this kind of 
collective excitation modes. We have shown that the 
calculated B(E2) and B(M1) in terms of the RPA cor- 
respond to those given by the conventional rotational 
model expressions, where the intrinsic quadrupole mo- 
ment and the effective ^-factors are calculated within the 
mean-field approximation. The link between the RPA 
and the rotational model expressions is given in the adi- 
abatic limit, where the precession phonon energy goes to 
zero. Then the rotor gn-1 actor is not a common factor 
any more, but depends on the configurations, especially 
on the occupation of the high-j decoupled proton orbital. 
Since the RPA formalism includes this effect properly, the 
calculated B(M1) values reproduces the experimentally 
deduced ones rather well. It is, however, noticed that 
the adiabatic approximation is not necessarily a good ap- 
proximation because the precession energies are not very 
small; more crucially, if a high-j decoupled orbital with 
Q = 1/2 is occupied, the approximation breaks down 
completely. Therefore, it is an important future task to 
examine how the non adiabatic effect plays a role in the 
realistic cases. More experimental data, especially B(E2) 
and B(M1) values, are necessary for this purpose. 

Finally, it is worth mentioning the similar RPA calcu- 
lations for the wobbling motion in the Lu and Hf region. 
We have presented the result in the recent papers pdlll2fl . 



Although we obtained the RPA solutions, which have ex- 
pected properties of the wobbling motion, the calculated 
out-of-band over in-band B(E2) ratios were smaller than 
the measured ones by about a factor two to three; this 
was the most serious problem in our microscopic calcu- 
lation. The measured ratio is almost reproduced by the 
simple rotor model. Both the out-of-band and in-band 
B(E2), which are vertical and horizontal transition dis- 
cussed in Sec. |n] are expressed in terms of the intrinsic 
quadrupole moments, Q20 an d Q22 .30] (or e.g. defor- 
mation parameters (£2,7)), combined with the wobbling 
phonon amplitudes. In the RPA wobbling formalism, on 
the other hand, the in-band transition is calculated by 
the intrinsic moments, while the out-of-band transition 
by the RPA phonon transition in Eq. (|23|l . Thus the un- 
derestimation of the B(E2) ratio above means that the 
RPA phonon transition amplitudes is smaller by about 
50—70% than the expected ones. 

The adiabatic approximation can also be considered 
for the case of the wobbling phonon 0] . Similar corre- 
spondence between the intrinsic moments and the RPA 
transition amplitudes, like Qo ~ (Qo)rpa in the present 
paper, is obtained with a non-trivial modification: There 
are two amplitudes related to the operators Q21 ^ an d 
Q 2 2^ m Eq. (|27|l . and the B(E2) is calculated by a linear 
combination of them with coefficients involving the three 
moments of inertia. Therefore, incorrect coefficients of 
amplitudes would make B(E2) values to deviate consid- 
erably, even though the adiabatic approximation is ap- 
plicable and two amplitudes are obtained in a good ap- 
proximation. There is, of course, another possibility that 
the adiabatic approximation itself is no longer valid. It 
should be noted that the wobbling excitation energies ob- 
served in Lu isotopes are about 200—500 keV, which are 
not small if translated to the transition phonon energy 
in the laboratory frame, w wo b + w rot ; see Eq. (|38ll . In 
the light of the present investigation, it may be possible 
that the RPA approach yields the correct magnitude of 
out-of-band transitions also for the case of the wobbling 
mode, because it actually does in the case of the preces- 
sion phonon bands. Thus, it is a very important future 
issue to examine whether or not the RPA wobbling for- 
malism can describe the observed B(E2) ratio in the Lu 
and Hf region. 
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